Occurrence of normal and anomalous diffusion in polygonal billiard channels 
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From extensive numerical simulations, we find that periodic polygonal billiard channels with angles which 
are irrational multiples of n generically exhibit normal diffusion (linear growth of the mean squared displace- 
ment) when they have a finite horizon, i.e. when no particle can travel arbitrarily far without colliding. For 
the infinite horizon case we present numerical tests showing that the mean squared displacement instead grows 
asymptotically as ?logf. When the unit cell contains accessible parallel scatterers, however, we always find 
anomalous super-diffusion, i.e. power-law growth with an exponent larger than 1. This behavior cannot be ac- 
counted for quantitatively by a simple continuous-time random walk model. Instead, we argue that anomalous 
diffusion correlates with the existence of families of propagating periodic orbits. Finally we show that when 
a configuration with parallel scatterers is approached there is a crossover from normal to anomalous diffusion, 
with the diffusion coefficient exhibiting a power-law divergence. 

PACS numbers: 05.45.Pq, 05.60.Cd, 02.50.-r, 05.40.Fb 
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I. INTRODUCTION 

Billiard models, in which point particles in free motion un- 
dergo elastic collisions with fixed obstacles, have simple mi- 
croscopic dynamics but strong 'statistical' behavior, such as 
diffusion, at the macroscopic level. The statistical properties 
of scattering billiards such as the Lorentz gas, which have 
smooth, convex scatterers, are well-understood: see e.g. id 
for a collection of reviews. Recently, however, models with 
polygonal scatterers have attracted attention, since they are 
only weakly chaotic, with for example all Lyapunov expo- 
nents being 0, but they still exhibit surprising statistical prop- 
erties in numerical experiments 

The dynamics in polygonal billiards depends strongly on 
the angles of the billiard table. Some rigorous results are 
available for angles which are all rational, i.e. rational mul- 
tiples of n: see |5] and references therein. In particular, an 
initial condition with a given angle cannot explore the whole 
phase space, since the possible angles obtained in the result- 
ing trajectory are restricted, so that the model does not have 
good ergodic properties. Previous numerical work has shown 
that if one angle is rational and the others irrational, then the 
results are sensitive to the value of the rational angle [3l|4j]. 

In this work we restrict our attention to the case in which 
all angles are irrational multiples of n. There are few rigorous 
results available in this case, but there is numerical evidence 
that the dynamics is mixing in irrational triangles |6]. 

We have performed extensive simulations of the transport 
properties of quasi-one-dimensional polygonal channel bil- 
liards; here we report detailed results for two classes of model. 
In both we find that transport generically corresponds to what 
is usually referred to as normal diffusion, i.e. linear growth 
of the mean squared displacement, when the particles cannot 
travel arbitrarily long distances without collisions (finite hori- 
zon). By 'generic' we mean that the set of models which do 
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not exhibit normal diffusion is 'small', for example having 
measure zero in the space of geometrical parameters. Sim- 
ilarly, for the infinite horizon case we find that transport is 
generically marginally super-diffusive, i.e. there is a logarith- 
mic correction to the linear growth law of the mean squared 
displacement. 



However, we also find exceptions to the generic behaviors 
described above. These exceptions occur when there are par- 
allel scatterers in the unit cell which are accessible from one 
another (see Sec. [HI] for the definition). In such cases we 
always find super-diffusion, i.e. the mean squared displace- 
ment grows like a power law with exponent greater than 1. 
Further, we demonstrate that as a parallel configuration is ap- 
proached, the transition from normal to anomalous diffusion 
occurs through a power-law divergence of the diffusion coef- 
ficient. We also show that the exponents we obtain are not 
consistent with a simple continuous-time random walk model 
of anomalous diffusion. 



Our results are consistent with previous observations of ir- 
rational models exhibiting normal 01 B 01 an d anomalous 
IE EE 03 Hill diffusion. Also, very recently, numerical re- 
sults corresponding to one of the models studied in this work 
were reported in 1 12], where it was also observed that parallel 
scatterers result in anomalous diffusion; the results we obtain 
complement and extend those of that reference. 



Energy transport in polygonal billiards placed between two 
heat reservoirs at different temperatures has also been a focus 
of attention OH 01- This process is closer to "color diffu- 
sion" than to heat conduction due to the absence of particle 
interactions and, hence, of local thermodynamic equilibrium 
fPjll . In these systems the properties of the energy transport 
process are known to be closely related to those of the particle 
diffusion process fl4lfl5ll . and for this reason we only study 
the latter. 
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FIG. 1: Unit cells of (a) the polygonal Lorentz model, and (b) the 
zigzag model. 



A. Models 

The unit cells of our models are shown in Fig.[Q the com- 
plete channel consists of an infinite horizontal periodic repeti- 
tion of such cells. The polygonal Lorentz model is a polygonal 
channel version of the Lorentz gas channel studied in 11611 . 
consisting of a quasi- ID channel with an extra square scat- 
terer at the center of each unit cell. By quasi-lD we mean 
that particles are confined to a strip infinitely extended in the 
x-direction but of bounded height in the y-direction. 

The polygonal Lorentz model can be simplified by elim- 
inating the central scatterer and flipping the bottom line of 
scatterers to point upwards, resulting in the quasi- ID zigzag 
model. Both models are designed to permit a finite horizon by 
blocking all infinite corridors in the structure, so that there is 
an upper bound on the distance a particle can travel without 
colliding with a scatterer; in scattering billiards this condi- 
tion was shown to be necessary for normal diffusion to occur 
fl7lfT3l . The zigzag model has also been studied in 1 12]. 

We remark that the zigzag model with parallel sides can be 
reduced to a parallelogram with irrational angles. This simple 
billiard geometry seems to have been neglected previously, 
although it is close to a model considered in 1 19] which was 
reported to exhibit a very slow exploration of phase space. 
We have found similar behavior in the parallel zigzag model 
(see also I12I0 . but it appears to be a peculiarity of this model 
which is not necessarily shared by other polygonal billiards 
exhibiting anomalous diffusion. 



B. Simulations 

In our simulations we distribute N = IQ 1 initial conditions 
(unless otherwise stated) uniformly with respect to Liouville 
measure, i.e. with positions distributed uniformly in the avail- 
able space in one unit cell, and velocities with unit speed and 
uniformly distributed angles. We consider only irrational an- 
gles for the billiard walls in each model, although we vary 
explicitly scatterer heights rather than angles, and to fix the 
length scale we take d = 1 (see Fig.Q . 

We study statistical properties of the particle positions x(t) 



as a function of time t, denoting averages at time t over all 
initial conditions by (•),. Of particular interest is the mean 
squared displacement (A* 2 ),, where Ax, :— x(t) — x(0), which 
is frequently used to characterize the transport properties of 
such systems. 



II. NORMAL AND MARGINALLY ANOMALOUS 
DIFFUSION 

In this paper, by normal diffusion we mean asymptotic lin- 
ear growth of the mean squared displacement (Ax 2 ), as a func- 
tion of time t , which we denote by (Ax 2 ), ~ t as t — > °°. This is 
characteristic of systems commonly thought of as 'diffusive', 
for example random walkers which take uncorrelated steps of 
finite mean squared length. 

Strongly chaotic dynamical systems, such as the Lorentz 
gas, are known to exhibit normal diffusive behavior [jl ; this 
arises due to a fast decay of correlations, even between neigh- 
boring trajectories, so that after a short correlation time the 
dynamics looks 'random'. It is rather surprising, then, that 
normal diffusion also seems to appear in non-chaotic polygo- 
nal billiards |3|,0,|7|] m which correlations do not decay expo- 
nentially. In these systems, the only source of 'randomization' 
arises due to divergence of neighboring trajectories colliding 
on different sides of a corner of the billiard: see e.g. |20[ for a 
recent attempt to characterise this randomization effect. 

In the following, we present sensitive tests which provide 
strong confirmation that, apart from exceptional cases, the 
decorrelation induced by the randomization mechanism men- 
tioned above is indeed sufficient to give rise to normal diffu- 
sion in polygonal billiards. In particular our tests distinguish 
the logarithmic correction (marginally anomalous diffusion) 
which arises generically in these systems when they have in- 
finite horizon, as it does in fully chaotic billiards. 

We will see that exceptions to these behaviors occur when 
the billiards present accesible parallel walls. In such cases we 
argue that long time correlations persist, related to the exis- 
tence of families of travelling periodic orbits, which always 
give rise to anomalous diffusion. 



A. Growth of second moment 

The growth of the mean squared displacement (Ax 2 ), as a 
function of time t is the most commonly used indicator of 
transport properties. In Fig. [2 we plot (Ax 2 ), for several non- 
parallel zigzag models, varying h\. For values of h\ below 
/i3 = 0.45 there is a finite horizon and we find that a linear 
regime is rapidly attained, suggesting normal diffusion. For 
values of h\ above this value, there is an infinite horizon; in 
this case, the mean squared displacement is still almost linear, 
although a small amount of curvature is visible. 

Normal diffusion corresponds to an asymptotic slope 1 of 
(Ax 2 ), on a double logarithmic plot. Figure [3] shows such a 
plot for the same data as in Fig. [2] and we see that indeed the 
slope is 1 for finite horizon and close to 1 for infinite horizon 
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FIG. 2: (Color online) (Ax 2 )/ as a function of t in the zigzag model 
with I12 = 0.77 and I13 = 0.45. There is a finite horizon only when 
hi < I13 = 0.45; curvature is barely visible in the infinite horizon 
cases. Statistical errors are smaller than the line width. 
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FIG. 3: (Color online) (Ax 2 ), as a function of t on a double logarith- 
mic plot in the zigzag model with the same parameters as in Fig. [2] 
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FIG. 4: (Color online) {Ax 2 ) t /t as a function of logf in the zigzag 
model with the same parameters as in Fig. [2] 



where z '■— Inf. Asymptotic linear growth of (Ax 2 ) t /t as a 
function of lnf thus indicates marginal anomalous diffusion, 
whilst asymptotic flatness corresponds to normal diffusion. 
A similar method was used in 12311 for a discrete time sys- 
tem, and we have checked that it also correctly finds marginal 
anomalous diffusion in the infinite-horizon Lorentz gas with 
continous-time dynamics. Figure|4]plots (Ax 2 ) t /t, confirming 
the normal/marginal distinction in finite and infinite horizon 
regimes. 

We can also fit (Ax 2 ), to the expected form at\nt + bt + 
c with parameters a, b and c, for example using a nonlinear 
least-squares method. For h\ = 0.71 we find a = 0.051; this 
agrees very well with the asymptotic slope 0.049 in Fig. 0] 
as it should do by equation d2. 1 1 . For h\ = 0.3 we instead 
obtain a = —0.0003, i.e. essentially zero (note that it cannot 
be negative), confirming normal diffusion. 



(e.g. slope 1.06 for h = 0.71). For models with irrational an- 
gles and no accessible parallel sides (see Sec.lllU. we always 
find a growth exponent close to 1, so that normal diffusion is 
generic in polygonal billiard channels. 

The above test does not, however, take account of logarith- 
mic corrections such as those found in the Lorentz gas with 
infinite horizon, where (Ax 2 ) r ~ t lnf ifTsIl : such a correction 
will appear on the above plot as a small change in the ex- 
ponent. A heuristic argument of l2lll leads us to expect that 
in any channel with an infinite horizon, including polygonal 
ones, we will have (at least) this kind of 'marginal' anomalous 
diffusion. We now present numerical tests which distinguish 
the logarithmic correction, showing that the diffusion is nor- 
mal in the finite horizon case and marginally anomalous in the 
infinite horizon case. 

In marginal anomalous diffusion we expect that (Ax ) t ~ 
at In f + bt + c for some constants a, b and c 12211 , so that 

(Ax 2 ), 

^-^-~az + b + ce- z , (2.1) 



B. Integrated velocity autocorrelation function 

For a system with normal diffusion we define the diffusion 
coefficient D as the asymptotic growth rate lim,^<>o ^(Ax 2 ), 
of the mean squared displacement. The Green-Kubo formula 
D = lim,^oo /oC(r)dT I24J1 relates D to the velocity autocor- 
relation function C(t ) = (v, vo): D exists, and the diffusion is 
normal, only if C(f ) decays faster than 1/f as f — * °°. C(f ) was 
studied in \% for a different channel geometry, but this func- 
tion is very noisy and it is not possible to determine its decay 
rate. 

Instead, following 1 25J we consider the integrated velocity 
autocorrelation function 

t 

R(t):= fc(T)dT = (v Ax t ), (2.2) 


where Ax r := x(f) — x(0) is the particle displacement at time 
f. This function is much smoother than C(f), and satisfies 
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FIG. 5: (Color online) Integrated velocity autocorrelation function 
R(t) as a function of log t in the (non-parallel) zigzag model with the 
same parameters as in Fig.[2|(/i2 = 0.77 and hj = 0.45). 



R(t) — > D as f — ► °° if the diffusion coefficient exists, whilst 
R(t) — lnf if C(t) decays as f -1 . 

Figure[5]plots R(t) as a function of lnf for several finite and 
infinite horizon non-parallel zigzag models. The flatness of 
R(t) in the finite horizon models indicates that the limit exists, 
and hence that the diffusion is normal, and contrasts with the 
asymptotic linear growth in the infinite horizon models, pro- 
viding further evidence of the marginally anomalous behavior 
in the latter case. 



C. Higher-order moments 



by 



We define the growth exponent y q of the gth moment ( \x\ q ) t 



t^°° logf 



(2.3) 



which ignores corrections to power-law growth and is a con- 
vex function of q; further, y q < q for all q, since the particles 
have finite velocity |26J. The exponents y q , and in particu- 
lar the diffusion exponent a = Y2, are measured using a fit to 
the long-time region of a double logarithmic plot of the rele- 
vant moment, assuming that the asymptotic regime has been 
reached. 

We find that the behavior of moments of high order (q > 4) 
is dominated by the extreme particles, i.e. by those which 
have the largest value of \x(t)\, giving data which is not repro- 
ducible between different runs. For this reason we eliminate 
the 5 most extreme particles from the average before calcu- 
lating the growth rate, resulting in data which is now repro- 
ducible. A similar procedure was used in 12/Tl . 

Figure |6] shows the exponents y q for the zigzag model with 
finite and infinite horizon. The low-order exponents satisfy 
Yq = l/2, and in both cases there is a crossover to a second lin- 
ear regime with slope 1, which occurs at approximately q = 3 
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FIG. 6: (Color online) Growth rate y q of qth moment (|x| 9 )f for 
finite horizon {h\ = 0.4) and infinite horizon (hi = 0.5) in the zigzag 
model with hi = 0.77 and A3 = 0.45 for 10 6 initial conditions. 



in the finite horizon case and q = 2 for infinite horizon; the lat- 
ter agrees with the result found in 1 26] for the infinite horizon 
Lorentz gas using a method which also applies here. Higher- 
order moments thus provide another method to distinguish be- 
tween normal and marginal anomalous diffusion in polygonal 
billiards. 

The observed qualitative change in behavior of the mo- 
ments corresponds to a change in relative importance between 
diffusive and ballistic effects 1 26, |28l |29[, and explains the 
previous observation that the Burnett coefficient, defined as 



the growth rate of the fourth cumulant (x 4 ) r — 3 (x 2 ) ( z , diverges 
as a function of time in polygonal billiards (Hy)]. 

The combination of the above tests provides strong evi- 
dence of the distinction between normal diffusion in finite 
horizon models and marginally anomalous diffusion in infi- 
nite horizon models. 



III. ANOMALOUS DIFFUSION 

Whilst we generically find normal diffusion, for certain ex- 
ceptional geometrical configurations this does not occur. The 
main exception is when there are accessible parallel scatter- 
ers in the unit cell of a model. In this case we always find 
anomalous (super-)diffusion, i.e. (Ax 2 ) t ~ t a with a > 1: see 
Fig- El Since the particles have finite speed, we always have 
a < 2, with the value 2 corresponding to ballistic motion. 
Sub-diffusion (a < 1) was reported in J^l in a model with 
one rational angle, but we are not aware of any model with 
irrational angles exhibiting sub-diffusion. 

We say that a model has 'accessible parallel scatterers' 
when there are trajectories which can propagate arbitrarily far 
by colliding only with scatterers which are parallel to some 
other scattering element in the unit cell. 

As an example, the results we report for the polygonal 
Lorentz model were obtained with a square central obsta- 
cle, that is, a scatterer whose opposite sides are parallel. 
This alone does not, however, imply that the system exhibits 
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FIG. 7: (Color online) Mean squared displacement (Ax 2 )? for par- 
allel models versus ( on a double logarithmic plot. The upper and 
lower lines are proportional to ; and f 2 for comparison. The remain- 
ing curves show, from top to bottom, the parallel zigzag model with 
h\ =0.1 and h-x, = 5.2, 10.2 (both shifted vertically for clarity), and 
the parallel polygonal Lorentz model with hi = 0.1, w = 0.15 and 
/j3 = 5.2, 10.2. Statistical errors are slightly greater than the width of 
the lines. 



anomalous diffusion, since it is not possible for a particle to 
propagate by colliding only with (copies of) this square, i.e. 
the scatterers making up the square are not 'accessible' from 
one another. 

Any particle must in fact also collide with the top and bot- 
tom scattering walls. Only if these walls are also accessibly 
parallel as defined above (i.e. have parallel partners in the unit 
cell) do we find anomalous diffusion. Furthermore, if we al- 
ter the central obstacle shape so that it no longer has parallel 
sides, then we still find anomalous diffusion provided that the 
scattering walls still have parallel partners which are accessi- 
ble from one cell to another. 

Figure[8] shows the behavior of the anomalous diffusion ex- 
ponent a as a function of I13 in the parallel zigzag and polyg- 
onal Lorentz models with finite horizon. The representative 
error bars shown give the standard deviation of the exponent 
over 100 independent simulations for the same value of /13, 
and indicate that the structure visible in the curve for the 
zigzag model is not an artifact. 

Figure |9l shows the growth exponent y q as a function of q 
in two models with anomalous diffusion. Again we find a 
crossover between two linear regimes, with y q of the form 1 28] 



Yq 



V\q, for q<q c 
q-V2, for q>q c . 



(3.1) 
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FIG. 8: (Color online) Diffusion exponent a as a function of A3 in 
the parallel zigzag model with h\ =0.1 and the parallel polygonal 
Lorentz model with h\ =0.1 and w = 0. 15, calculated using 10 5 ini- 
tial conditions evolved until the presumed asymptotic regime. (Note 
that at hj, = 1 = d there is a rational angle in both models.) Repre- 
sentative error bars show the standard deviation of the exponent over 
100 independent simulations. 
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FIG. 9: (Color online) Growth exponent y q of moments for the par- 
allel polygonal Lorentz and zigzag models with h\ = 0.51, /12 = /13 = 
10.0 and w = 0.21. 



A. Time evolution of densities 

A feature typical of normal diffusive systems is that the 
probability distribution of displacements (or equivalently of 
positions) converges in the long-time limit, when rescaled by 
y/i, to to a normal distribution 11611 . We find this behavior 
numerically in both of the classes of polygonal billiard chan- 
nels studied here, in the cases for which the mean squared 
displacement grows linearly, and we conjecture that this holds 
in general, providing a stronger sense in which polygonal bil- 
liards can be considered diffusive. 

For anomalous diffusion with (Ax-) t ~ t a (05 > 1), we no 
longer expect the central limit theorem to hold. The inset of 
Fig.llOlshows the probability density p t (x) of positions x(t) in 
a parallel polygonal Lorentz model. It exhibits a striking fine 
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FIG. 10: (Color online) Rescaled demodulated densities in the 
polygonal Lorentz model, compared to a Gaussian with the cor- 
responding variance. The inset shows the original and demodu- 
lated (heavy line) densities for t = 100. The model parameters are 
hi = 0.1, h 2 = h 3 = 0.45, w = 0.2, d=l. 
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FIG. 11: (Color online) Scatterplot representing the joint density 
yf{x,t) of laminar stretches in the zigzag model with hi = 0.1, ho, = 
0.3 and h2 = /?3 +Ah. The straight lines with slope 1 correspond to 
the maximum possible speed. 



structure, similar to that found in normally diffusive Lorentz 
gases and polygonal billiards in 0, where it was shown that 
the principal contribution to this fine structure arises from a 
geometrical effect, namely that the amount of vertical space 
available for particles to occupy varies along the channel. As- 
suming that the dynamics has good mixing properties, it was 
shown that defining the demodulated density f, (x) by 

f t (x):=p t (x)/h(x) (3.2) 

eliminates much of the fine structure. Here, h(x) is the height 
of the available space in the billiard channel at horizontal co- 
ordinate x, normalised by the area of the unit cell so that the 
integral J d d h(x)dx over the unit cell is 1. This demodulated 
density is also shown in the inset of Fig. ^| although the 
demodulation is not as effective as in the normally diffusive 
cases studied in [16], due to the weaker mixing properties of 
the current model. 

Rescaling the demodulated position density / ( (x) at time t 

by 

f t {x):=t a l\f t {xt a l 2 ) (3.3) 

gives a bounded variance as t — ► °° Hill . The main part of Fig- 
ure ^3 shows demodulated densities rescaled in this way for 
different times, compared to a Gaussian with the correspond- 
ing variance. The rescaled densities converge at long times to 
a non-Gaussian limiting shape, as has previously been found 
in other systems with anomalous diffusion ||30|] . 



B. Mechanism for anomalous diffusion 

1. Continuous -time random walk approach 

For the parallel zigzag model, Figure^] demonstrates the 
presence of long laminar stretches, i.e. periods of motion 



during which a particle's velocity component parallel to the 
channel wall does not change sign. Typical trajectories ex- 
hibit such behavior, sometimes over very long periods of time. 
Figure ^2 shows a scatterplot representing the joint distribu- 
tion of laminar stretches of duration t with horizontal 
displacement x along the channel. (To compute these data, 
the channel was divided into sections lying between adjacent 
maxima and minima in order to determine the wall direction). 

Thus, we may attempt to describe the dynamics of this 
system as a continuous-time random walk (CTRW), whose 
steps are the laminar stretches: see e.g. 13 ll 13211 . To this 
end, we note that Figuremshows that the distribution y/(x,f) 
is concentrated on diagonal lines emanating from the origin, 
which leads us to try an approximation of the form y/(x,f) = 
i 5(|jcj —vt) \jf(t) for some speed v. This version of the CTRW 
was termed the velocity model in l32ll . and describes motion 
at a constant velocity for a time t in the direction x; after 
each stretch the direction is randomised and a new step is 
taken. With this form for y/(x,t), the long-time growth of the 
mean squared displacement depends on the asymptotic decay 
rate of the marginal distribution \ff(t) of step times ll32ll : if 
\j/(t) ~ ?~'~ v as t — ► 00, then the variance <7 2 (f) behaves like 
<J 2 (t) ~ f 3 ~ v when 1 < v < 2, while we have normal diffusion 
a 2 (f) — t for v > 2. 

Figure [T2l shows the tail region of *P(f) := 1 — /q y/(T)dT. 
From the top curve we find a long-time power-law decay 
x P(f) ~ f~ L58 for a particular case having parallel scatterers. 
For close-to-parallel configurations, the tail of the distribu- 
tion follows that for the parallel case, but with an exponen- 
tial cut-off, as shown in the inset of Fig. ^] According to 
the CTRW model, this cut-off gives rise to normal diffusion 
in the asymptotic long-time regime. This change of behav- 
ior is in qualitative agreement with our observations. Unfor- 
tunately, the value of the tail exponent for the parallel case 
corresponds to v = 1.58 and hence to a mean squared dis- 
placement (7 2 (f) ~ f 3 ~ L58 — r 1 - 42 , which does not agree with 
the numerically observed growth rate <7 2 (t) ~ t . 
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FIG. 12: (Color online) Tail of ¥(?) for the parallel zigzag model. 
The top curve is for hi =0.1 and hj, = 0.3. The lower three curves 
are for h l = 0.49, h 3 = 5.1 and Ah = 0, 0.0001 and 0.001 (top to 
bottom). The inset is a semilogarithmic plot of the Ah = 0.001 case, 
showing exponential decay. 



Although the above results rest on a rather gross approxi- 
mation to the joint distribution \j/(x,t), we doubt that a more 
general CTRW model, incorporating information on the com- 
plete \j/(x,t), would substantially improve the agreement. The 
reason is that correlations between consecutive laminar trajec- 
tories have an important effect not accounted for by CTRW 
models. Indeed, for a given trajectory, we find that the lengths 
of consecutive long laminar stretches are often nearly equal, 
indicating that the trajectory repeats very closely its previous 
behavior many times. 



2. Propagating periodic orbits 

On the other hand, our results indicate that it is the presence 
of accessible parallel scatterers that gives rise to anomalous 
diffusion. We believe that propagating periodic orbits, i.e. or- 
bits which repeat themselves with a spatial displacement, may 
provide a mechanism for this behavior. Such orbits are more 
prevalent when there are parallel scatterers, since in this case 
it is much easier for a particle to regain its original angle of 
propagation after a given number of bounces, and, then, pos- 
sibly repeat its previous motion. In addition, periodic orbits in 
polygonal billiards occur in families 13311 : several such fami- 
lies are shown in Fig. 1131 

Furthermore, trajectories with initial conditions which are 
close in phase space to those of a propagating periodic orbit 
will shadow it, and hence propagate, for a long time, with a 
longer shadowing time for closer initial conditions. Anoma- 
lous diffusion should then result from a balance between the 
ballistic propagation of the periodic orbits themselves, the 
long-lasting ballistic motion of shadowing orbits, and the dif- 
fusive motion of other trajectories. 

We have, however, not yet been able to account analytically 
for anomalous diffusion in this way. For instance, the mere ex- 
istence of propagating orbits is not enough to give rise to the 



FIG. 1 3 : (Color online) Families of trapped and propagating periodic 
orbits in parallel systems. 

kind of anomalous behavior we observe. Indeed, families of 
propagating orbits exist in the corridors of the infinite-horizon 
periodic Lorentz gas, but these give rise only to marginally 
anomalous diffusion (Ax 2 ) ~ tint | 18] (see Sec.[n). Thus, al- 
though the set of propagating periodic orbits must have mea- 
sure zero in phase space, since otherwise their ballistic nature 
would result in overall ballistic transport with (Ax 2 ) ~ t 2 , they 
must be plentiful in some sense, in order to give rise to anoma- 
lous transport. For example, such orbits could be dense at least 
in some parts of phase space. 



IV. CROSSOVER FROM NORMAL TO ANOMALOUS 
DIFFUSION 

Since anomalous diffusion occurs only for special geomet- 
rical configurations, it is of interest to study the crossover from 
normal to anomalous diffusion as such configurations are ap- 
proached. We begin by considering the polygonal Lorentz 
channel, in which we fix all geometrical parameters except 
for hi, which we vary as h% = hi +Ah. 

Figure [H] shows a double logarithmic plot of the mean 
squared displacement in a polygonal Lorentz model for val- 
ues of Ah tending to 0, i.e. approaching the parallel config- 
uration; a linear fit to the bottom curve has been subtracted 
from each curve for clarity. As Ah tends to 0, the curve for the 
parallel configuration (Ah = 0) is followed for progressively 
longer times before a crossover occurs to asymptotic diffusive 
behaviour, which appears on the figure as a zero slope. Defin- 
ing the crossover time T c (Ah) as the intersection of the initial 
anomalous growth with a straight line fit to the asymptotic 
normal growth, we find that T c scales like T c ~ (A/;) -0,37 as 
Ah^O. 

The inset of Fig. [3] shows the diffusion coefficient D as 
a function of Ah, obtained from the asymptotic slope of the 
mean squared displacement. Although Fig.ll4lshows that the 
asymptotic linear regime has not been reached for the small- 
est values of Ah, thereby underestimating the diffusion coeffi- 
cient, we obtain a straight line on a double logarithmic plot of 
D versus Ah, giving the power-law behavior 

D{Ah) ~ |A/!p 015 as |A/i| 0, (4.1) 
with the same exponent for both positive and negative values 
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FIG. 14: (Color online) Double logarithmic plot of {Ax 2 ) t as a 
function of time for the polygonal Lorentz model with hi = 0.1, 
/?3 = 0.45, w = 0.2 and I12 = hj, + Ah. Deviations from an asymp- 
totic linear fit to the lowest curve are shown for Ah = (uppermost) 
and Ah = 10~ m , m = 2, . . . , 11 (bottom to top). The inset shows a 
double logarithmic plot of the diffusion coefficient D as a function of 
Ah for the polygonal Lorentz model with hi = 0. 1 for I13 = 0.3, 0.45, 
and for /i 3 = 0.45 with Ah < 0. 



1() 4 




10-' 10° 10 1 10 2 10 3 



u = \Ah\ v t 

FIG. 15: (Color online) Data collapse of \Ah\ va (x 2 (Ah)), as a func- 
tion of u = |A/i| v f for the zigzag model with hi = 0.1, hi = 0.3, 
which has a ~ 1.81 (upper curve) and the polygonal Lorentz model 
with hi = 0.1, h3 = 0.45, with a = 1.40. In both models the 
value v = 0.37 was used, and curves are shown for Ah = 10 _m , 
m = 4, . , . , 11. The solid lines show the extreme behaviors of the 
scaling function 14.31 . 



of Ah and for two (close) values of ho,. For this model, the 
diffusion exponent in the parallel (Ah = 0) case is a ~ 1 .40. 

A similar crossover to anomalous diffusion occurs near to 
the parallel configuration in the zigzag model (not shown), 
and we conjecture that such crossover behavior is generally 
found in polygonal billiards near to configurations exhibiting 
anomalous diffusion. The rate of growth of the diffusion coef- 
ficient for Ah 7^ depends, however, on the diffusion exponent 
for Ah — 0: for the zigzag model with hi =0.1 and /13 = 0.3, 
we find D ~ |AA| -0 ' 2 , with a diffusion exponent a ~ 1.81 
for the parallel case; both of these exponents differ from the 
results in the polygonal Lorentz model. 

The growth of the diffusion coefficient close to a parallel 
case thus depends on the anomalous diffusion exponent found 
in that parallel case. Indeed, the relation between the ex- 
ponents characterizing the crossover from normal to anoma- 
lous transport can be obtained from the following simple 
scaling argument. Define the crossover time exponent V by 
T c ~ \Ah\- y . We have 



(x\Ah))t 



jD(Ah)t, for t>T c 
\t a , 



forf < T c , 



(4.2) 



so that continuity at T c gives D{Ah) ~ 7;"- 1 ~ \Ah\ v( - l - a \ 
which is in good quantitative agreement with our results. 

Furthermore, this argument implies that plots of 
\Ah\ va {x 2 (Ah)) t as a function of 11 := \Ah\ v t should col- 
lapse onto a scaling function of the form 



0(h) 



for u > 1 
for «Cl. 



(4.3) 



The data collapse for both the zigzag and the polygonal model 
near a configuration of parallel scatterers is shown in Fig. 1151 



Interestingly, while the transport exponent a is sensitive to 
the geometrical details of the models, the best collapse was 
obtained for the same value of v = 0.37 in both cases. 

Finally, it should be noted that the above calculations all 
refer to irrational angles; the situation when the value of Ah 
renders an angle rational is unclear. 



V. DISCUSSION AND CONCLUSIONS 

We have demonstrated that the type of diffusive behavior 
exhibited by periodic polygonal billiard channels depends on 
geometrical features of the unit cell. Based on our results, we 
conjecture that sufficient conditions for a periodic polygonal 
billiard channel to exhibit normal diffusion are the following: 
(i) all vertex angles are irrationally related to n\ (ii) the billiard 
has a finite horizon; and (iii) there are no accessible parallel 
scatterers in the unit cell. 

When there is an infinite horizon, diffusive transport is re- 
placed by marginal anomalous diffusion with (Ax 2 ) ~ t lnf, 
while if there are accessible parallel scatterers, then anoma- 
lous diffusion with {Ax 2 } ~ t a , a > 1 is observed. 

In the zigzag model we also have evidence of anomalous 
diffusion (e.g. with exponent 1.16 for one particular model) 
with irrational angles chosen such that one angle is twice the 
other, i.e. with a rational relation between them, in disagree- 
ment with the result found in fl2ll (where a rather small num- 
ber of initial conditions was used). For other angle ratios we 
do not have conclusive data, but the exponents are seemingly 
close to 1, while for the polygonal Lorentz model we find nor- 
mal diffusion in similar situations. For the zigzag model we 
must thus also exclude the possibility of rationally-related an- 
gles from the above conditions for normal diffusion. 

We remark that in |7|], propagating periodic orbits were 
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conjectured to be related to anomalous diffusion in a system 
with one rational angle, but no reason was given for their ex- 
istence. Actually, it is possible to include that system into our 
picture by unfolding the rational angle, which gives rise to an 
equivalent system with parallel scatterers. 

We believe that the explanation of anomalous diffusion 
should be found in terms of such propagating periodic orbits, 
which are much more prevalent in the presence of parallel 
scatterers. We further showed that there is a crossover from 
normal to anomalous diffusion as a parallel configuration is 
approached, with the diffusion coefficient having a power-law 
divergence. We hope to achieve a quantitative description of 
both of these points in the future. 
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